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ABSTRACT 

In the context of a cosmological study of the bulk flows in the Universe, we present a detailed study of 
the statistical properties of the kinematic Sunyaev-Zel'dovich (kSZ) effect. We first compute analytically the 
correlation function and the power spectrum of the projected peculiar velocities of galaxy clusters. By taking 
into account the spatial clustering properties of these sources, we perform a line-of-sight computation of the 
all-sky kSZ power spectrum and find that at large angular scales (/ < 10), the local bulk flow should leave a 
visible signature above the Poisson-like fluctuations dominant at smaller scales, while the coupling of density 
and velocity fluctuations should give much smaller contribution. We conduct an analysis of the prospects 
of future high resolution CMB experiments (such as ACT and SPT) to detect the kSZ signal and to extract 
cosmological information and dark energy constraints from it. We present two complementary methods, one 
suitable for "deep and narrow" surveys such as ACT and one suitable for "wide and shallow" surveys such 
as SPT. Both methods can constraint the equation of state of dark energy w to about 5-10% when applied to 
forthcoming and future surveys, and probe w in complementary redshift ranges, which could shed some light 
on its time evolution. This is mainly due to the high sensitivity of the peculiar velocity field to the onset of 
the late acceleration of the Universe. We stress that this determination of w does not rely on the knowledge of 
cluster masses, although it relies on cluster redshifts and makes minimal assumptions on cluster physics. 

Subject headings: cosmic microwave background - large scale structure of universe 
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1. INTRODUCTION 

The new generation of ground-based high-resolution cos- 
mic microwave background (CMB) ex periments ( e. g., th e 
Atacama Co s mology Telescope 3 [ACTKosowskv (2003); 
iFowler et al J J2005h1 and the South Pole Telescope 4 
[SPTRuhl et al. (2004)]), are designed to scan with very high 
sensitivity and arcminute resolution the microwave sky. Their 
main goal i s the study of the ther mal Sunyaev-Zel'dovich 
(tSZ) effect (S unvaev & Zerdovichll980l) : the change of fre- 
quency of CMB photons due to inverse Compton scattering 
by hot electrons. Such hot electron plasma are known to 
be found in clusters of galaxies, and should also be present 
in larger structures, such as filaments and superclusters of 
galaxies. This scattering translates into a redshift indepen- 
dent distortion of the CMB black body spectrum, making the 
tSZ effect an ideal tool to probe the baryon distribution in the 
large scales of our Universe at different cosmic epochs. How- 
ever, this is not the only effect of an electron plasma on the 
CMB radiation. If a cloud of electrons is moving with some 
bulk velocity with respect to the CMB frame, then Thomson 
scattering by these electrons will imprint new (Doppler in- 
duced) temperature fluctuations on the CMB photons. This 
is the so called kinematic S unyaev-Zel'dovich effect (kSZ, 
ISunvaev & Zeldovichl ( 119721) '). which is spectrally indistin- 
guishable from the intrinsic CMB temperature fluctuations. 

Although the kSZ effect is typically an order of magnitude 
smaller than the tSZ in clusters of galaxies (and for this rea- 
son much harder to detect), it encodes precious cosmological 
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information since it depends on the peculiar velocity field. 
Indeed, kSZ measurements can yield valuable information 
about the large-scale velocity field, the evolution of the dark 
matter potential, and the growth of fluctuations. The study 
of large scale velocity fields (cosmic flows) has been an ac- 
tive research area in the nineties. There have been numerous 
attempts to measure bulk flows using the large-scale distri- 
bution of galaxies and their peculiar velocities, and to place 
constraints on the matter power spect rum or the Univ erse 
matte r density, (see revie ws of Strauss & Willick ( 119951) and 
ICourteau & DekeJ (1200 1) and references therein). However, 
it became clear that these measurements had to be corrected 
for systematic errors, such as the biases introduced when cal- 
ibrating the distances of the galaxies under study, or the non 
linear components of the velocities of those objects. With kSZ 
observations, by using clusters as tracers of the velocity fields, 
one is more confident to probe larger (less non-linear) scales. 

While there have been no kSZ detections to date, upper 
limits on the pe culiar velocities of i ndividual clusters have 
been reported by Benson et alj ( 120031) . Such a difficult mea- 
surement could in principle be hampered by other effects 
such as non-line arities and the co mplicate d physics of the 
intra-cluster gas. Nagai, Kravtsov, & Kosowskv (2003) have 
shown that the kSZ is not diluted by the internal v elocity 
dispersion in the intracluster gas. Ma & Fry (2002) calcu- 
lated the temperature fluctuations produced b y the kSZ in 
the no n-linear regime using the halo model. iBenson et alJ 
( 120031) showed that the signal-to-noise of a kSZ measure- 
ment should be distance independent and sugg ested com- 
bining signals from different redshifts. Holder (2002) and 
Agha nim et alJ ll2004l) ha ve discussed h o w to e xtract the kSZ 
signal from maps, while Schafer et al. (2004) used N-body 
simulations to build template s of kSZ maps in the context 
of the Planck mission, and Kashlinsk v & Atrio-Bara ndela 
(2000) study the possibility of extracting the kSZ dipole from 
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CMB surveys covering a large fraction of the sky. Due 
to its weak signal, indirect detection of the kSZ has been 
proposed throu gh cross-correlation techn iques with weak- 
lensing j Pore, Hennawi, & Spergell 120041) or old galaxies 
llDeDeo. Sp ergel & Trac 2005). As we shall see below, a mea- 
surement of the kSZ effect would be very valuable, since it 
would not only allow us to measure bulk flows of clusters of 
galaxies and test the predictions of the standard model, but 
also provide additional constraints on cosmological parame- 
ters, especially on the equation of state of dark energy. 

In this paper we compute the correlation function and the 
power spectrum of the kSZ effect. This requires modeling of 
the peculiar velocity field and the cluster population. We ex- 
plore the prospects for future CMB experiments to measure 
the kSZ correlation function, and find that ACT-like experi- 
ments should be able to detect the kSZ-induced CMB vari- 
ance at high (~ 12cr) significance level. Further, we study 
the dependence of the kSZ correlation function on the cosmo- 
logical parameters, and show that it can be used to measure 
the equation of state of dark energy (w) if the redshifts of the 
clusters of galaxies detected in CMB surveys are available. 
We find that, with the SALT follow up of ACT data or with 
a wider but shallower (SPT-like) survey, the w parameter can 
be constrained with an accuracy of 8% for an ACT scan of 
400 square degrees. This error should scale inversely with the 
square root of the covered area and hence becomes 5% for a 
1,000 square degree area. These determinations of w do not 
rely on the knowledge of cluster masses, but do rely on clus- 
ters redshifts and reasonable assumptions on cluster physics. 

Unless otherwise stated, throughout this work we shall 
asume a LCDM cosmological model ( ISoergel et al.l l200"3h 
with n m = 0.3, fl A = 0.7, h = 0.72, and cr 8 = 0.88. The paper 
is organised as follows: in Section 2 we study the statistical 
properties of the projected peculiar velocity field and provide 
an analytical expression for the correlation function of 
projected velocities. In Section 3 we compare the kSZ and 
tSZ effects, and discuss the strategy to enhance the probabil- 
ities of detecting the former. In Section 4 we compute the 
correlation function and the power spectrum of the kSZ, both 
when we consider only a given set of clusters present in a 
survey or the whole celestial sphere. In Section 5 we outline 
two methods to estimate the kSZ effect in future CMB cluster 
surveys and in Section 6 we explore the dependence of kSZ 
measurements on cosmological parameters making particular 
emphasis on w. We conclude in Section 7. 



2. THE CORRELATION FUNCTION OF LINE-OF-SIGHT 
LINEAR PECULIAR VELOCITIES 

While the measurem ent of the kSZ effect of individual 
cluster s is difficult (e.g.. Aghanim et al. (2001); Benso net al.l 
(2003)), in this paper we address the prospects for statistical 
detection of cluster peculiar velocities in future CMB surveys. 
This requires the knowledge of the ensemble properties of 
the velocity field traced by the galaxy cluster population, to 
which we devote the current section. For clarity and future 
reference, a statistical description of the linear velocity field 
and related quantities is given in Appendix A. 

Throughout this paper, we shall assume that the measured 
velocity field obeys linear theory. However, as noted by 
Colbe rget al.l (120001) . this is not completely fulfilled by 
cluster velocities, since clusters are peculiar tracers of the 
large scale matter distribution and show biased velocities 



compared to t he exp ectations provided by the linear theory. 
Colbe rget al.l (12000) found that this bias was typically a 30% 
- 40% effect. Although several attempts have been made 
to model this boost in terms of the un derlying density field 
(Sheth & Diaferio 2001; Hamana et al. 2003), in subsequent 
sections it will be accounted for by simpl y increasing the 
cluste r velocities by a factor b v =1.3 dSheth & Diaferiol 
2001). The goal of this paper is to present a theoretical 
calculation of the detectability of the signal and the fore- 
casted signal-to-noise: detailed comparison with numerical 
simulations are left to future work. Full treatment can be 
implemented only with the help of numerical simulations 
matched to a given observing programCsee lPeell (120051) for a 
recent study). We must stress that, although some modelling 
of non-linear effects must be included in this study, clusters 
are the largest virialised structures known in the Universe, 
and probe much bigger scales than galaxies. Therefore we 
must expect them to be significantly better tracers of the 
linear velocity field. Furthermore, as we shall see below, the 
peculiar velocity estimator (the kSZ effect) does not depend 
of distance, which avoids the need to use redshift-independent 
distance indicators, as opposed to other peculiar velocity 
estimators. 

On large, linear scales, the density and peculiar veloci- 
ties are related through the continuity equation: dS^/dt = 
-/k • \it I a, where a and k are the scale factor and comoving 
Fourier mode, respectively. The peculiar velocity of a cluster, 
as probed by its kSZ effect, can be interpreted as the linear 
peculiar velocity field smoothed on comoving scale R which 
corresponds to the cluster's mass M via 
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where p is the background matter density. The kSZ effect 
is sensitive to the line-of-sight component of the velocity, 
but under the assumption that the velocity field is Gaussian 
and isotropic (which should be satisfied in the linear regime) 
the three spatial components of the velocity field must be 
statistically independent. Moreover, the power spectrum 
must completely determine the statistical properties of the 
velocity field. Thus, in a given cosmological model, the linear 
velocity field power spectrum (which in turn is related to 
the matter power spectrum) should univocally determine the 
angular correlation function (and angular power spectrum) of 
the line-of-sight cluster velocities. 

In linear theory, the velocity dispersion smoothed over 
spheres of comoving radius R (corresponding to a given clus- 
ter mass M) is given by 

2 



al,(R,z)=[H(z) 



dV s 



dz 



dk k z 



{k)l W(kR)\\ 



2n 2 k 2 



(2) 



where W(kR) is the Fourier transform of the top hat win- 
dow function, H(z) is the Hubble function, T>${z) is the linear 
growth factor and P,„(k) is the present day linear matter power 
spectrum, (the power spectra at any redshift will be denoted 
as P$s(z,k) = T>j(z) P m (k)). Hence the power spectrum of the 
velocities is: 
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where T> v = H{z)dT)$ / dz is the velocity growth factor. When 
computing T>$ for different Dark Energy models, we used the 
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FIG. 1 . — Peculiar velocity growth factor for three different cosmologies: a 
ACDM model (fi t = 0.7, Q m = 0.3, thick solid line), a fi m = 0.3, flat w = -0.6 
model (dotted line), and Q m = —1 /3 and w = — 1/3, (dashed line). 



analytical fit provided bv lLinderl (120051) : 



g(a) = exp / dloga 
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where = T>s{a) ja gives the deviation of the growth factor 
from that of a critical (f2„, = 1) universe, and 7 is given by 

7 = 0.55+&[l+w(z=l)], (5) 

with b = 0.05 if w > -1 and b = 0.02 otherwise. 

The growth of the velocity perturbations with redshift may 
provide useful cosmological constrai nts such as constraints on 
the e quation of state of dark energy ( DeDeo, Sperg eT& Trad 
2005). This is illustrated in Fig.Q, where we show the 
redshift evolution of the velocity growth factor for three 
different cosmological models: a ACDM model (Oa = 0.7, 
Oi = 0.3, thick solid line), a flat universe with il m = 0.3 and 
dark energy equation of state parameter w = -0.6 (dotted 
line), and another flat model with Q m = 0.3 and w = -1/3, 
(dashed line). Due to the k 2 factor in the denominator of 
eq-O, the signal is weighted by the largest scales, making 
this probe relatively insensitive to the smoothing scale and 
therefore to clusters mass. Indeed, if the dependence of a vv 
versus mass is approximated by a power law, then one finds 
that for the concordance model <j vv oc M~ 013 . 

Having this in mind, we compute here the angular corre- 
lation fu nction of th e line-of-sight (LOS) cluster velocities. 
Note that lPeell (12005) takes a different approach to this calcu- 
lation. Assuming that we can measure the LOS component of 
the peculiar velocity of a cluster, we compute the quantity 



C w (0 12 ) = ( vCxO-m v(x 2 )-n 2 



(6) 



where ni and n 2 denote two different directions in the sky, 
"connecting" the observer to the cluster positions X; and x 2 
and 6 \2 denotes the angle between ni and n 2 . We refer the 
reader to Appendix B for the detailed derivation and here we 
report the final expression for this correlation function: 
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In this equation, the factor T\ is given by 
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and X\,X2 are the (comoving) distances to the clusters, 
(without loss of generality we have used the convention that 
x\ > x 2 ). Here ji(x) denote the spherical Bessel functions and 
the summation must take place only over even values of Z; R t 
and R2 refer to the linear scales corresponding to the masses 
of each cluster. Note that we recover eq.(|2} in the limit of 
On — > and x\ — > x%. 

Fig.(|2^) shows the behaviour of C vv vs d\2 in the concor- 
dance ACDM model for a couple of 10 14 M clusters when 
they are both placed at z=0.005 (solid line), both placed 
at z=0.1 both (dashed line) and when once cluster is at 
z=0.1 and the other at z=l (dotted line). In the first (clearly 
unrealistic) case, the clusters are so close to us that both are 
comoving in the same bulk flow with respect to the CMB 
frame, giving rise to the dipolar pattern shown by the solid 
line. In the case where both clusters are at z ~ 0.1 (dashed 
line), their correlation properties are strongly dependent on 
9\2, since this angle defines the distance between the clusters. 
For small angular separation, the clusters are still relatively 
nearby, and hence their peculiar velocities are correlated, but 
this correlation dies as the separation of the clusters increases. 
The angular distance at which the correlation drops to half its 
value at zero separation is ~ 10°, corresponding to roughly 
40 h~ l Mpc. Finally, if clusters are very far apart from each 
other (thick dotted line), their peculiar velocities are not 
correlated. 

An alternative way to present the correlation properties of 
the projected peculiar velocities of the clusters is the veloc- 
ity field angular power spectrum. In Appendix B, we invert 
C vv {9\2) into its angular power spectrum C"', i.e., 



Q,(0u) = J2 ^-CrP,(cos9 l2 ). 

1 

We find that 

CT=^-j-(/B/-i + (/+l)S, + i), 
where the Bi coefficients are defined by 
k 2 dk 
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Fig. 0b displays the power spectra for two cases: two very 
nearby clusters (both at z=0.005, solid line), showing an 
almost dipolar pattern, and two relatively far away clusters 
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FIG. 2. — Correlation function (top) and corresponding power spectrum 
(bottom) of the projected peculiar proper velocities. The solid lines observe 
the case where the two clusters are placed at z = 0.005 from the observer. 
Dashed lines show the case in which clusters are further, z = 0.01. Finally, 
the dotted line considers the case where one cluster is located at z = 0. 1 and 
the other cluster is much further (z= 1). Note the lack of correlation in this 
case. 



(both at z=0.1, dashed line), case in which the power is 
transferred to higher multipoles. The power spectrum for the 
clusters placed at z=0.005 and z=0. 1 is zero. 



ative velocity to the CMB frame. An observer will only be 
sensitive to the radial component of the electron peculiar ve- 
locity, so the expression for the change in brightness temper- 
ature becomes (Sunyaev & Zel'dovich 1968) 



ST 



kSZ 
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= dl n e (l) <jj 



(12) 



Here, we have assumed that the peculiar velocity is the 
same for all electrons, r stands for the optical depth, and 
n is a unitary vector giving the direction of observation. 
This process must take place in two contexts: (1) When 
the intergalactic medium becomes ionized by the high 
energy photons emitted by the first stars, inhomogeneities 
in the electron velocity and density distributions generate a 
kSZ signal which is known as the Ostriker-Vishniac effect. 
Despite of the large size of ionized structures encountered by 
the CMB photons, the electron density contrast is relatively 
small, and further, we do not know where in CMB maps to 
look for this signal because we do not know the location of 
the ionized bubbles which formed during reionization. The 
amplitudes and angular scales at which this effect should 
be visible are model dependent, but can be as high as a 
few m icroK in the multipole range I > 2000, (Sant os et alJ 
120031) . (2) Clusters of galaxies at lower redshift leave a more 
easily detectable signal. Their high electron density can give 
rise to values of t as high as 10~ 3 - 10~ 2 , and their peculiar 
velocities should be close to 300 kms" 1 at z = 0, which 
together can produce temperature fluctuations of the order of 
1-10 /iK. The clusters high optical depth is also responsible 
for large spectral distortions of the CMB generated through 
the thermal Sunyaev-Zel'dovich (tSZ) effect. The tSZ effect 
is typically an order of magnitude larger than the kSZ and 
introduces frequency dependent brightness temperature 
fluctuations which, in the not relativistic limit, change sign at 
218 GHz. Therefore by combining observations in bands at 
frequencies lower and higher than this cross frequency, it is 
possible to obtain the cluster position and to characterize the 
tSZ cluster signal. Once the cluster position is identified their 
kSZ contrib ution should be accessible at 218 GH z. However, 
as noted b y ISeheal. Kosowskv & HolderJ (12005 ) and earlier 
bv lHolderi 12001 7 even with measurements in three different 
frequencies it may not be possible to obtain a clean estimate 
of the kSZ effect for a single cluster. We shall show below 
that we are not interested in a very accurate kSZ estimate 
for a given cluster, but on unbiased estimates on our entire 
cluster sample. 

In the next subsection we make a detailed comparison of 
the amplitude of the kSZ and the tSZ effects. 



Here we have characterised the projected peculiar velocity 
field at cluster scales. Next, we address the study of the kSZ 
effect and its comparison with the tSZ effect. 

3. THE SUNYAEV-ZEL'DOVICH EFFECTS 

3.1. The kinematic Sunyaev-Zel 'dovich Effect 

The Kinematic Sunyaev-Zel'dovich effect describes the 
Doppler kick that CMB photons experience when they en- 
counter a moving cloud of electrons. Since this is simple 
Thomson scattering, there is no change of the photon fre- 
quency, and hence it leaves no spectral signatures in the CMB 
blackbody spectrum. Therefore, this effect is solely deter- 
mined by the number density of free electrons and their rel- 



3.2. Comparison of the kSZ and the tSZ effects in clusters of 

galaxies 

In what follows we shall describe the galaxy clus- 
ter population by ado p ting the model presented in 
IVerde. Haiman & Spergel (2002). This model is based 
upon the spherical collapse description of galaxy clusters 
dGunn & Gotll972i) . and assumes that clusters are isothermal 
and their gas acquires the virial temperature of the halo. The 
halo mass and redshift distribution is approximated by the 
formalism presented in Sheth & Tormer] (11999ft . We refer to 
Verde. Haiman & Spergel (2002) for further details in this 
modelling. 
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Both kSZ and tSZ effects can be written as integrals of 
some function K.(r) along the line of sight crossing the cluster, 
weighted by the electron density: 

— =g(x) J dla T n e (l)JC{l) (13) 

For the non-relativistic tSZ g(x) = (x cothx/2 - 4), 
x = hv/ksTo is the adimensional frequency in terms of 
the CMB monopole To, K. = kgT? / (m e c 2 ) where n e , T e , m e are 
the electron density, temperature and mass respectively, or 
the Thomson cross-section and ks the Boltzmann constant. 
For the kSZ, g(x) = 1, K = — n- v/c, where n is a unitary 
vector pointing along the line of sight, and v is the cluster 
peculiar velocity. If clusters are perfectly virialised objects, 
then one must expect the scaling T e oc M 2 / 3 . However, 
following IVerde. Ha iman & Srjereel ( 2002), the cluster 
model should leave some room for some deviations from 
such scaling, which can be due to internal (non-linear) cluster 
physics, deviations from purely gravitational equilibrium, 
preheating, etc. Our parametrization adopts T e oc M 1 '^, where 
£ = 1 .5 for a perfectly virialized cluster, but spans the range 
from 1.5 to 2 in the literature. We shall present results for 
£ = 1 .5, but the sensitivity of our results on £ is minimal. 
At the same time, the radial peculiar velocity dispersion of 
the clusters, (y/(—n- v/c) 2 ) = a vv decreases very slowly with 
mass (~ M~ 013 ), so, at fixed redshift we must expect the ratio 
kSZ/tSZ to be bigger for low mass clusters, which are more 
numerous. Regarding the redshift dependence, clusters tend 
to be denser at earlier epochs (the product n e r v scales roughly 
as (1 + z) 2 ), so we must expect larger tSZ and kSZ amplitudes 
at higher redshifts. However, the dependence of the function 
JC is different in each case: while clusters are hotter at earlier 
times (T e oc (1 + z)), velocities have not had so much time to 
grow as at present epochs (<r vv oc 1 /y/l+z), and hence the 
ratio kSZ/tSZ (oc 1/(1 + z) 3/2 ) decreases with redshift. 

We show these scalings explicitely in Fig.(|3}, where the thin 
lines evaluate the tSZ and the kSZ at z = and thick lines 
correspond to z = 1 . The solid lines provide the amplitude 
of the cluster-induced tSZ fluctuations at 222 GHz, which, as 
we shall see below, will be taken as our effective frequency 
after accounting for the tSZ relativistic corrections and the 
effects related to the finite spectral width of the detectors. The 
dashed lines provide the expected rms amplitude of the kSZ 
effect in clusters: we see that, at z = 0, clusters of masses 
below ~ 10 15 Mq should produce more kSZ than tSZ flux at 
222 GHz, at least by a factor of a few. This mass threshold 
decreases at z = 1, but so does the typical halo mass at such 
redshift (clusters of ~ 10 14 M© are very rare objects), so again 
we must expect the kSZ to dominate over the tSZ. 

4. THE POWER SPECTRUM OF THE KSZ EFFECT 

4.1. Model of the cluster population 

We next study the kSZ signal generated by the entire pop- 
ulation of clusters of galaxies by computing its two second- 
order momenta, i.e., the correlation function and the power 
spectrum. For this, it is first necessary to have a model to 
describe the population of galaxy clusters in the Universe. 
We shall adopt the hierarchical scenario in which small scale 
overdensities in the Universe become non-linear and collapse 
first, and then merge and give rise to bigger non-linear struc- 
tures. The abundance of haloes of a given mass at a given 
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FIG. 3. — Comparison of the tSZ (solid lines) and kSZ (dashed line) effects 
as a function of clusters mass. Thick lines refer to z=l (where most of the 
clusters are located for ACT-like surveys), whereas thin lines correspond to 
z=0. For the majority of clusters, the kSZ flux will be a few times bigger than 
the tSZ flux at 222 GHz. 



cosmic epoch or redshift is given by the clusters mass func- 
tion. We will adopt the Sheth & Tormen (hereafter ST) mass 
function ( Sheth & Tormer] ri999l) . which will be denoted by 
h(M,z) = dN(M,z)/dV(z) and provides the average number 
density of haloes of mass contained between M and M+dM 
at redshift z. It must not be confused with n(M,x), which is 
the actual number of haloes in that mass range at position x. 
The latter can be understood as a random variable, the for- 
mer as its mean. As it will be useful later, we first compute 
the mean number of haloes present in two volume elements 
centered at X] and X2: 

(n(Mi ! Xi)«(M 2 ,x 2 ))=«(Mi,zi)n(M 2 ,Z2) + 

6 D (M, -M 2 )5 3 D (xi -x 2 )n(M u zi) + 

(A[n(M 1 ,x 1 )]A[n(M 2 ,x 2 )]). (14) 

In this equation, zi and zi are the redshift corresponding to po- 
sitions Xj and X2 respectively. The first term in the right hand 
side of the equation is merely a constant, but will have its rel- 
evance, because it will couple with the velocity field as we 
shall see below. The next term containing the Dirac deltas ac- 
counts for the (assumed) Poissonian statistics ruling the (dis- 
crete) number density of sources, and will be referred to as 
the Poissonian term. In the third term, A[n(Mi,xi)] stands 
for the deviation with respect to the average halo number den- 
sity due to the environment, i.e., due to large scale overden- 
sities, which condition the halo clustering. Therefore, this 
third term describes the spatial clustering of haloes, which 
is a biased tracer of the spatial clustering of matter. In the 
extended Press-Schecter approach it can be shown that the 
power spectrum and the correlation function of haloes and un- 
derlying matter are merely proportional to each other over a 
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wide rang e of scales. This is commonly expressed by a bias 
factor <Mo & White! 19961) . so that 



&*(r) = b 2 (M,z) W) ; Phhik) = b 2 (M,z)P m (k). 



(15) 



Here £,hh,£,m and Phh(k),P m (k) stand for halo-halo and 
matter correlation function/power spectrum, respectively. 
Therefore, the third term in the RHS of eq. (!14> equals 

n(M u zi)n(.M2,Z2MM 1 ,ziMM2,Z2)^ m (xi-x 2 ). 

A parallel approach consists in writing the halo mass func- 
tion as a function of some linear scale matter overdensity 
6 = (p—p)/p. The number of haloes atx is then approximated 

as 



n(M,x) = n(M,z) + r]+ 



dn(M,z) 



OS 



S(x) + 0[S 2 ], (16) 



,5=0 



where r\ is a random variable which introduces the Poisso- 
nian behaviour of the source counts. In the extended Press- 
Schechter formalism, it turns out that dn(M,z)/dS\ g _ Q coin- 
cides with the bias factor b(M,z), and by using this it is possi- 
ble to reproduce eq. H4i from ea.( ll6> . This justifies neglect- 
ing all higher-order powers of 8 in eq. dl6l >. We have taken 
S to be in linear regime, but, since £/,/,(>) ~ b 2 ^ mm {r) down to 
scales comparable to the halo size, we shall use this formalism 
down to halo scales. 

4.2. A Line of Sight approach for the kSZ effect 

We next write the temperature anisotropics induced by the 
kSZ in clusters of galaxies as an integral along the line of 
sight: 



ATI 



kSZ 



To 



[n]: 



dr Tj 



W gas 
j 



(r-rj). (17) 



Here, r is the comoving radius integrated to the last scatter- 
ing surface, and the sum over the index j represents a sum on 
all clusters; fj denotes the opacity in the center of the cluster 
(f = a<jjn e . c with n e>c the central electron number density and 
a the scale factor), and the window function Wj ;fls (r-r ; ) de- 
notes the gas profile of the cluster. Alth ough we should adop t 
some realistic shape for this profile ( Komatsu & Seliak 2001 ), 
we have adopted a simple Gaussian window with scale radius 
equal to the virial radius of the cluster. This is justified since, 
in Fourier domain, at scales bigger than the cluster size, the 
window function is merely equal to the volume occupied by 
the gas, regardless of the shape of the gas profile. This step 
simplifies our computations significantly, and does not com- 
promise the accuracy in the relatively big scales (bulk flow 
scales) in which we are interested. Note that even if the sum 
is made over all clusters, only the clusters being intersected 
by the line of sight will contribute to the integral. This sum 
can be re-written first as an integral and then as a convolution, 



A7i 



kSZ 



To 



[n]: 



dr / dM x 



v(M,y)-n 



n[M,y] W gas [M 7 r-y] 



dyf(M,z) 



(f(M,z) [-**£^ ) „[M])*W« fli [M,z] 



dr / dM x 



M. (18) 



Note that central optical depth and the window function 
have an intrinsic dependence on redshift. This is due to the 
fact that clusters formed at higher redshift tend to be more 
concentrated. The symbol * denotes here convolution in real 
space. 



4.3. The all sky correlation function and power spectrum of 
the kSZ effect 

If two different lines of sight are now combined to estimate 
the angular correlation function, then one obtains that 



7b T 
n(M u ri)[ v(M 1 ,r,)-n 1 )n(M 2 ,r 2 )( v(M 2 ,r 2 )-n 2 



(19) 

Next we plug ea.(ll6> and make use of the Cumulant Ex- 
pansion Theorem. We must note as well that Poissonian fluc- 
tuations will be assumed to be independent of 5, and that for 
Gaussian statistics, the 3- and 4-point functions are zero, just 
as (v) and (5) . Therefore we are left with a sum of products 
of 2-point functions of the form 



n(Mi,Zi)<4 v (M h zi) + 
n(M u zi)n(M 2 ,z 2 )( ( v(Afi,ri)-ni ) ( v(M 2 ,r 2 )-n 2 



dn{Mi,zi) 



OS 



dn(M 2 ,z 2 ) 



;5=0 



OS 



i=0 



(<Kri)5(r 2 ))( v(Afi,r!)-ni v(M 2 ,r 2 )-n 2 ) + 



<5(ri) v(M 2 ,r 2 )-n 2 )( S(r 2 ) v(Mi,ri)-n 



S(r{) [\(Mi,r{)-ni 



<5(r 2 ) v(M 2 ,r 2 )-n 2 



(20) 



Note that the last term is a constant. We refer again to 
Appendix A where the cross terms (S(\ ■ n)) are studied. In 
Appendix C we provide a explicit computation of the power 
spectra arising from each of the terms considered in eq.( !20t . 
Since the last term introduces no anisotropy, the kSZ power 
spectrum is the sum of four contributions: a Poisson term 
(Poisson), a term proportional to the velocity-velocity cor- 
relation (vv term), a term proportional to the product of the 
velocity-velocity correlation and the density -density correla- 
tion (vv - dd term), and a term proportional to the density- 
velocity squared (dv-vd term). Therefore, 

c, = cf + c; T + cf d ~ vv + cf v ~ vd . (2i) 

Fig.Q displays each of the terms in ea. (l21> : the thick solid 
line is the Poisson term, whereas the thick dashed line corre- 
sponds to the vv term. Although the latter decreases rapidly 
with increasing I, at the large scales it dominates over all other 
terms, reflecting the presence of the local bulk flow. In an at- 
tempt to simplify the expressions for these two terms given in 
Appendix C, we have found the approximate integral 



dz 



dV(z) 
dz 



dMn(M,z)(T 2 JM,z) 



\yt(M,z)\ 



(22) 
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FIG. 4. — Different components of the all-sky kSZ power spectra: the 
thick solid line shows the Poisson term, which dominates over the other terms 
except in the low / limit. In this large scale range, the vv generated by the local 
bulk flow is the one introducing most power, (dashed line). Filled circles and 
diamonds show semi-an alytical approximations for these two Poisson and 
vv, respectively (see eas. 1221231 ). Of lower amplitude, the thick dotted line 
shows the sum of the dd-vv and the dv-vd terms, (thin dotted lines, note that 
the latter is negative). 




Angle (0 / degrees) 

FIG. 5. — Cluster-cluster kSZ correlation functions for a sample of clusters 
with masses above 2 X 1O 14 M0. When forming cluster pairs we require pair 
constituents to be at similar redshifts. If we consider all redshifts simultane- 
ously, we obtain the solid line. The dashed lines correspond to differential-vn- 
redshift cluster-cluster correlation functions: the thick dashed line considers 
only clusters around z ~ 0.01, whereas the intermediate dashed line corre- 
sponds to z ~ 0. 1 and the thin dashed line to z ~ 1 . 



for the Poisson term, and 

cp 



1 



dM—n(M,z) yi(M,z) 

4"7T 1 



(23) 



for the vv term. y/(M,z) is the Fourier transform of the cluster 
profile in the sphere, 



2tt0 v exp-/(/+ 1)0;/2, 



(24) 



and P vv (k, z) is the Ic and z dependent velocity power spectrum. 
6 V is the angular size of the cluster virial radius. Filled circles 
for the Poisson term and diamonds for the vv term provide a 
comparison of these approximations with the exact integrals. 
Although the amplitudes and slopes are not too disimilar, 
the approximated vv term is particularly innacurate at large 
angular scales (I < 10). The approximation for the Poissonian 
term seems to predict a somewhat correct amplitude at low 
/, but a shallower slope, which translates into a smaller 
amplitude at high multipoles. 

The dd-vv and dv-vd terms are shown by the thin dotted 
lines. We must remark that the dv-vd term is negative, and we 
are plotting its absolute value. The sum of both is given by the 
thick dotted line. The sum of these two terms is particularly 
hard to detect, since it never dominates, not at large scales (it 
is about a factor of 20 below the vv term), nor at small scales, 
where it is well below the Poisson term. 



4.4. The cluster-cluster kSZ correlation function 

Future high resolution multi frequency CMB experiments 
like ACT or SPT can provide kSZ estimates on those regions 
of the sky where clusters of galaxies have been identified via 
tSZ. Therefore one could attempt to measure the velocities 
correlation function by computing the kSZ correlation func- 
tion in this restricted set of pixels: 

C c k ! iz cl (9) = J2Y, T(Mi < > */> , z j2 ) x 



(25) 



v(M, 1 ,z ;i )-n 1 v(M, 2 ,z >2 )-n 2 



Since the signal comes from clusters which velocities are 
correlated, we should consider pairs close in redshift. The 
solid line in Fig.(|5j shows the correlation function to be mea- 
sured from all galaxy clusters from z = upto z = 4 more mas- 
sive than 2 x l0 H h~ 1 M Q , according to the standard ACDM 
cosmology and the Sheth-Tormen mass function. Members 
of the cluster pairs must be within Az = 0.01. Since the kSZ 
amplitude per cluster is typically of few tens of micro Kelvin, 
the zero-lag correlation function can be as high as ~ 150 (/iK 
) 2 , and drops to one half of this value at 9 ~ 2° - 3°. The thick 
dashed line shows the correlation function for clusters located 
at z — 0.01. As the redshift increases (z — 0.1 medium thick- 
ness dashed line; z — 1, thin dashed line), the amplitude at 
zero lag increases (clusters are more concentrated) and the 
correlation angle decreases. The solid line shows the redshift- 
integrated cluster-cluster correlation function: the signal is 
dominated by high redshift clusters, more concentrated and 
more numerous per unit solid angle. 

Diaferio, Sunv aev & Nusserl J2000I) pointed out a poten- 
tially important non-linear aspect that is not included in our 
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modeling: in very massive superclusters, the kSZ effect shows 
typically a dipolar pattern, plausably caused by the encounter 
of two oppossite bulk flows at their common attractor's po- 
sition. Since our model predicts no dipolar pattern at scales 
of a few degrees, such scenario is not accounted for by our 
approach. Hence, in a realistic application the core of such 
overdense regions should be excised from the analyses. Such 
massive structures however are very rare and form at very late 
epochs, and their exclusion should not compromise the anal- 
ysis presented here. 

5. CAN THE KSZ EFFECT BE MEASURED? 

In this section we outline two different procedures to 
extract the kSZ signal from future high-resolution and high- 
sensitivity CMB experiments. The procedures presented here 
may be sub-optimal but our aim is to quantify the relative 
importance of difference sources of error and to roughly 
forecast the expected signal to noise for those experiments. 
We defer the developement of an optimal procedure to future 
work. We shall try to extract the kSZ sign al in a s tatistical 
sense: while previous works (e.g., lA ghanim et alJ |2004)) 
have addressed the difficulty of separating the kSZ effect 
from potential contaminants (tSZ, radio-source emission, 
infrared galaxies, CMB, etc) in a given cluster, our approach 
will consist on combining the signal coming from subsets 
of clusters in such a way that the contribution of the noise 
sources averages out. As we shall see, the success of this 
procedure will rely on the precision to which the average 
properties of the potential contaminants are known. 

The first approach, which we shall refer to as method (a), is 
based on measurements of the kSZ flux and its redshift evo- 
lution. Its sensitivity to cosmology increases with redshift, 
and, for this reason it is suited for "deep and narrow" survey 
strategies. Here we will use ACT's specifications. Our second 
approach (method (b)) uses the ratio of kSZ and tSZ induced 
temperature anisotropies. This ratio is particularly sensitive to 
w at z^0.8 and since cosmic variance for the peculiar velocity 
field is more important at low redshifts, this method is more 
suitable for "wide and shallow" survey strategies. Here we 
will use a survey with specifications similar to those of ACT, 
but covering 4,000 square degrees and thus with higher noise. 
These specifications are not too dissimilar from those of SPT, 
assuming that accurate photometric redshifts can be obtained 
for all SPT clusters. 

5.1. Probing the kSZflux at high redshift 

In what follows, we shall use specifications for ACT one 
year data (Gaussian beam with FWHM equal to 2 arcmins, 
noise amplitude lower than 2 fiK per beam, and a clean 
scanned area of 400 square degrees). ACT observes in three 
bands: 145, 220 and 250 GHz. We will concentrate on the 
220 GHz channel, although some knowledge will be assumed 
to be inferred from the other bands. For instance, the 250 
GHz channel will be very useful to estimate the level of 
infrared-galaxy emission at lower frequencies. Likewise, the 
145 GHz channel should be critical when characterizing the 
tSZ flux from each cluster. For simplicity, we will concentrate 
on the variance of the kSZ signal, but it is easy to estimate the 
kSZ angular correlation function introduced in subsection 4.4. 
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FIG. 6. — Typical amplitude of CMB residuals remaining after attempting 
to approximate the CMB average within a circular patch of angular radius 6 r 
(given in abscissas) by the CMB average computed within a ring surrounding 
the patch of width 10% the patch radius. For most of clusters located at 
z > 0.3, these residuals will typically be of 3 - 4 fiK amplitude. 



SALT 5 . Cluster detection via tSZ effect will provide targets 
for optical observations. Alternatively, optical cluster identi- 
fication should be possible up to z ~ 1 from SALT'S multi- 
band imaging us i ng algorithms such as those developed in 
iKim et al.N2002l) : iMiller et alJ J2005I) . Hence, a direct com- 
parison can be made with tSZ detected cluster sample. SALT 
spectroscopic follow up will enable to obtain the cluster red- 
shifts. 

The method is as follows: 

• For every detected galaxy cluster, we take a patch of ra- 
dius equal to one projected cluster virial radius. We can 
use the edge of significant tSZ emission at other fre- 
quencies to define this radius. We compute the mean 
temperature within this patch, and draw a ring sur- 
rounding it, of width, say, 10% of the virial radius. 
We next compute the mean temperature within this ring 
and subtract it from the mean of the patch. This oper- 
ation should remove most of the CMB contribution to 
the average temperature in the patch, but will unavoid- 
ably leave some residuals, which will be denoted here 
as 5T™ b . These residuals will have two different con- 
tributions: the first coming from the inaccurate CMB 
subtraction, the second one being due to instrumental 
noise residuals, 

((^cm B ) 2 ) = (^L) + -^. (26) 

['beams 

Nbeams is the number of beam sizes present in the ring, 
and N 2 is the instrumental noise variance. If we denote 
the patch of radius the cluster virial radius as region 1, 
and by region 2 the ring surrounding it, it can be easily 



The strip covered by ACT should also be surveyed by 



5 SALT'S URL site: |http : / /www . salt . ac . za/| 



9 



proved that the first term in the RHS of the last equation 
reads 



1 



subs I 



AOi J AO 



dnidn 2 C(ni,n 2 ) + 



dn 1 dn 2 C(ni,n 2 )- 



asi 2 J An 



A^iA0 2 



dn l dn 2 C(n u n 2 ). (27) 



ASli J An 



Afii and AS1 2 denote the solid angles of the patch 
and the ring, respectively, whereas ni and n 2 denote 
directions in the sky and C(ni,n 2 ) is the CMB angular 
correlation function evaluated at the angle separating 
the directions ni and n 2 . The rms fluctuations intro- 
duced by this residual are plotted in Fig.(|6j: although 
it can be as high as a few tens of microK for nearby 
clusters subtending 30 - 40 arcmins, their effect 
reduces to a 3 - 4 /iK for clusters of a few arcmins 
size, which correspond to most of clusters at z > 0.3 
to be detected by ACT-like CMB experiments. This 
contaminant will be dominant over the instrumental 
noise contribution. 



• We assume that, by observing the tSZ amplitudes at 145 
GHz and 250 GHz, it is possible to provide an unbiased 
estimate of the mean tSZ contribution to the cluster 
patch at 220 GHz. After subtracting the tSZ and CMB 
components estimates, the temperature in the patch cor- 
responding to cluster i can be written as 



* = ^ B : + 8T tSZ , i+ N + ST'f ZJ+ T ksZti . 



(28) 



5T,sz,i denotes the tSZ residuals after substracting the 
estimated tSZ amplitude. N accounts for contribution 
of instrumental noise to the patch average, and ST^' Z , 
for the residual contribution of in ternal velocities. 
According to Diaferio et al.l J2005I) . we assume a 
typical rms for internal velocities of one third the bulk 
flow velocity expected for each cluster, and random 
orientation (sign); thus we assume ST^' zj rms to be 
1 /3 of the bulk-flow induced kSZ amplitude. The first 
four quantities in the RHS of ea. (l28> have zero mean, 
hence <5,- is an unbiased estimate of the kSZ amplitude 
in cluster i. 



following kSZ variance 6 estimator at redshift zi along 
the direction n: 



P*k /(") = E ^7 $>A)tfA) y 



i<] 



KJ 



(29) 



ilj stands for the j-th cluster's solid angle. We remark 
that the kSZ flux depends solely on the cluster's mass 
and peculiar velocity, and so do our angle-integrated 
kSZ temperature anisotropics. If we denote by Ni and 
Nj the number of cluster members in bins I and J, re- 
spectively, then the number of cluster pairs that can be 
formed by combining these two mass bins is given by 
Nij; N u = NiNj if / 4 J and N,.j = N,(N, - l)/2 for the 
same bin (/ = J). The indexes i and j run for individ- 
ual cluster members in each mass bin. wjj is a weight 
factor, which we define as 



N,Nj 



'i u j 



N, 



mST^ B ) 2 +(ST lS z) 2 + ((6T>" s ' z )2) +N 2 ]; 
Nj 



(30) 



with the subscripts / and J evaluating the brackets in 
the corresponding mass bins. The estimator of eq.j29> 
provides a weighted measurement of the kSZ variance, 



kSZ. I 



(n)) = 



E/<y w i,J (Fksz, ijFksz, u) 



(31) 



with Fksz, 1 j the expected kSZ amplitude at redshift zi 
and mass equal to that corresponding to bin /. Its formal 
error is given by 

A 2 [F 2 SZJ (n)] = 



2x 



S/<y W U (Fksz, i.iFksz, u) 



E 



i<j 



1J.L 



WLjr(FkSZ. IjFkSZ. ll) 



wij,M(Fksz, ijFksz, i,m) 



1J.M 



KJ 



KJ 



• We now combine estimates of the kSZ coming from 
different clusters in order to estimate the kSZ cluster- 
cluster correlation function, (Tksz.iTkszj) ■ As men- 
tioned in a previous Section, we are interested in com- 
puting this correlation function by considering pairs of 
clusters of similar redshift. Note that we avoid squaring 
the kSZ estimate of the same cluster (i = j), since resid- 
uals would not average out and would introduce a bias 
in our kSZ variance estimates. Let us assume that the 
bulk flows occupy a typical scale (named here as co- 
herence scale) d co h on the sky. Within this patch along 
direction n, we can sort all clusters in mass and redshift 
bins. Let /, J, L and M denote the mass bins of clusters 
in a common redshift range centered at zi- By building 
cluster pairs upon all members of these bins, and com- 
bining different mass bins, it is possible to compute the 



' KJ 



(32) 



where W//^, = wuNl/ctI. Note that the first term in 
the RHS of this equation is not sensitive to the number 
of clusters within the coherence patch. Such term, 
containing the squared kSZ expectations for mass 
bins / and J, is associated to the (assumed) intrinsic 
Gaussian nature of the kSZ fluctuations, and it corre- 
sponds to the cosmic variance contribution. It will thus 
scale like the inverse of the survey area when different 

6 We integrate temperature fluctuations in the cluster's area, obtaining a 
quantity of units ^K strad. This integrated temperature fluctuations are es- 
sentially proportional to flux fluctuations, and for this reason they will be 
denoted by F^sz- One must keep in mind, however, that their units are not 
Janskys, but microKelvin X stereoradian. 
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coherence patches are combined in the analysis. The 
forth term is exclusively due to observational errors, 
whereas the second and third terms are hybrid: they 
show contributions from both the intrinsic uncertainty 
of the velocity field and observational errors. If Nb 
denotes the number of mass bins, and we take the 
weights, the cluster number and the F^sz equal for all 
mass bins (07 = a, Nj = N and Fksz.u = Fkszj for every 
mass bin /), then it can easily be proved that the forth 
term scales roughly as cr 4 /N /(Nb(Nb + 1)), that is, the 
squared variance expected for each cluster over the 
total number of pairs that can be formed. The second 
and third terms, in this limit, yield a 2 F^ sz t /(N(Nb + 1)), 
which scales inversely to the number of clusters. In 
our case, our kSZ flux estimates will be limited by the 
cosmic variance term. On the other hand, because of 
having very few objects, we considered one mass bin 
only. 



• Finally, we combine estimates from the N co h — 
^vfsky/Ocoh different projected coherent regions in the 
sky, where f s u v is the fraction of the sky covered by the 
CMB experiment, yielding 



TP l — 

r kSZ, 1 — 



E 



A 2 ,(n)L 



E 



1 



A 2 [^ z , ,00]" 



with an uncertainty 



A 2 t^k/(n)]' 



(33) 



(34) 



Note that the size of the coherent patch must depend of 
redshift: a bulk flow extenting up to 20 /r'Mpc at z=l 
subtends a degree on the sky, whereas if it is at z~ 0.05 
then it subtends around 8 degrees. 

We take the effective noise N to have a typical amplitude of 
5 /iK per beam 7 , and it accounts for both instrumental noise 
and the confusion noise associated to unresolved sources. 
STgjig is computed as stated in eq.j27l>. and contributes 
typically with a few microK. Regarding the tSZ residuals, 
ST,sz contains the contribution coming from relativistic 
tSZ corrections, and power leakage associated to the finite 
spectral width of the detectors. We conservatively assume 
that the amplitude of these residuals (remaining after the tSZ 
substraction) is typically the non-relativistic tSZ temperature 
increment expected at 222 GHz. We approximate the beam as 
a Gaussian of FWHM equal to 2 arcmins, and take 400 square 
degrees (f^ ~ 10~ 2 ) as the clean sky region covered by ACT. 
ACT's sensitivity limit is set to cluster above 2 x 10 14 /i _1 M Q , 
and the total number of clusters above this threshold predicted 
in this region of the sky by our Sheth-Tormen mass function 
is roughly 4,400. However, for this analysis we only use 
relatively bright and big clusters: the product of their angular 
size and tSZ temperature decrement at 145 GHz must be 
bigger than 160 /iK arcmin 2 , and this requirement decreases 
considerably the amount of available clusters. 

As a result of this section, in Fig.Q we show our expecta- 
tions of ACT's sensitivity on the kSZ variance when all clus- 
ters are grouped in the redshift bins z h , and = [0.02, 0.4, 0.8, 

7 The noise level for ACT is expected to be per beam. 
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FIG. 7. — Expected amplitude of the squared kSZ flux for ACT, after group- 
ing all cluster in six redshift bins. The solid line line corresponds to a ACDM. 
The dotted and the dashed lines display the kSZ variances for universes with 
Dark Energy components having w equal to -0.6 and -1/3, respectively 



1.2, 1.6, 2.0]. Solid, dotted and dashed lines correspond to a 
ACDM model(^ m = 0.3, Qa = 0.7), a flat universe with dark 
energy equation of state w = -0.6, and w = -1/3, respectively. 
We report error bars for the ACDM model: we drop the first 
point at z = 0.02 (which is dominated by cosmic variance) 
and focus in the high redshift range: the signal-to-noise ra- 
tio shows a maximum at z=0.8, and beyond this redshift the 
error bars start to increase due to the lack of massive clusters. 

5.2. The kSZ/tSZ ratio 1Z 

In this subsection, we address the study of the ratio of 
the kSZ and tSZ effects. When referring to temperature 
anisotropies, the kSZ and tSZ effects are integrals weighted 
by the cluster electron density along the cluster diameter, and 
are sensitive to the size and/or the concentration parameter 
of these objects. The ratio of the kSZ and the tSZ, however, 
should cancel these dependencies out to great extent, and 
provide a cleaner view of the cluster's temperature and 
peculiar velocity. 

Following our cluster model, we investigate the behaviour 
of the statistic 1Z = STusz /5T,sz in the cluster population. If 
the evolution of the cluster temperature is well described by 
the spherical collapse model, we shall find then that 7Z should 
be a good estimator of the dark energy equation of state w 
at recent epochs. But since cosmic variance becomes impor- 
tant at low redshifts (smaller volume for a given solid angle), 
in this case we shall compute our expectations for an experi- 
ment with sky coverage close to f s fy = 0.1 (4,000 square de- 
grees). To compensate for the wider survey area we assume 
a noise of 10 /iK per arcminute squared and that clusters up 
to z = 0.8 can be detected and resolved. We assume that the 
frequency coverage enables contaminant subtraction and ex- 
traction of clusters tSZ signal. We also assume that follow 
up observations will yield redshift for all observed clusters. 
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These specifications are not too dissimilar to those of the SPT 
telescope when combined with photometric follow up. 

We shall see that these are precisely the requirements 
needed to obtain cosmological information from 1Z. However, 
since in this case the kSZ signal is divided by the tSZ temper- 
ature decrement (measured at say, 145 GHz), we must be very 
careful with the noise contribution to the denominator of 1Z. 
Our approach will be to conservatively consider only clusters 
whose integrated tSZ temperature decrements are larger than 
160 /iK-arcmin 2 at 145 GHz. This implies that tSZ errors 
will be typically 5% - 10% of the estimated tSZ temperature 
decrement, and that they can be treated perturb atively as er- 
rors in the kSZ estimation in the numerator. Therefore, our 
model for the estimation of 1Z in a single cluster will be given 
by: 



STSmb i + ST tSZti +N+5T$ z ,+ e T kszJ + T m ,i 



(35) 



Most of the terms of this equation are defined exactly as in 
ea. J28> . The only new term is eTksz.i, which accounts for the 
extra error introduced by the uncertainty in the denominator, 
Ttsz.i- Here T t sz,i is the absolute value of the tSZ decrement of 
the cluster at 145 GHz and e was taken to be a normal random 
variable of rms 0.05 (5% error in T,szj which reflects into 5% 
error in Tisz.i)- We neglect the correlation of the errors in T,sz j 
with those of Tksz.i and take e independent of the other noise 
sources. 

We now write the analogue to eq.(!29> as 



(36) 



i<j ij ' r<j 
where the weights are defined as 

NiNj 



Wrj 



2„2 



a, a 



mc e M B ) 2 + msz) 2 +wi$ z ) 2 )+e 2 tf 



,sz) + <(^) 2 >++^ 2 



mzj + N 2 ]/ 



(37) 



The estimate of 72/ (n) in a redshift band zi and direction n, 
together with its uncertainty A 2 [72/(n)] can be obtained from 
eas. ( 1311321 by simply replacing F^ sz ,(n) by 72/(n). Similarly 

the express ions for 72/ and A 2 [72./] can be obtained from 
eqs. EED , 

In Fig.([8j we plot the ratio 1Z at different redshifts as it 
would be seen by an experiment like SPT. As before, we have 
assumed that only clusters above 2 x 10 14 /t'Mq can be seen, 
and grouped all clusters in the redshift bins z'i and = [0.02, 0.4, 
0.8, 1.2, 1.6, 2.0]. The lines refer to the same cosmologi- 
cal models as in Fig.Q. Note the similarity of this plot and 
Fig.Q. We see that 1Z is sensitive to cosmology at much 
lower redshifts (z< 0.8) than F^sz, and for the sensitivity anal- 
ysis following in the next Section we have only used the first 
three redshift bins. 



N 

\ '*.* 
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1 .0 

Redshift (z) 



2.0 



FIG. 8. — Expected kSZ/tSZ ratio TZ as measured by SPT for six redshift 
bins. The solid line corresponds to a ACDM model while the dotted and the 
dashed lines are for flat models with w equal to -0.6 and -1/3, respectively. 



6. KSZ SENSITIVITY TO MEASURING COSMOLOGICAL 
PARAMETERS 

We can now explore the dependence of the two statistics in- 
troduced in the previous Section on cosmological parameters. 



We define the x as 



E 



2; 



Isz, / Q\sz, I 



^IQlszj] 



(38) 



where Q 2 ksz , can either refer to the angle-averaged kSZ 
anisotropy (F) or the kSZ/tSZ ratio (72). The likelihood is 
thus C oc exp-|x 2 . We estimate errors using the Fisher ma- 
trix approach. In principle the parameters that enter in the 
analysis are J7 m , w, the fraction of cluster mass in the in- 
tra cluster medium /icm, the present-day normalization of the 
matter density fluctuations erg and the reduce Hubble constant 
h. In practice both methods are non sensitive to /icm, &8 and 
h separately for a fixed number of detected clusters, but on 
their combination in the form of an overall amplitude of the 
kSZ signal, which we shall asign to an amplitude parameter A 
alone. 

We assume that A is redshift-independent (or that its 
scaling in redshift can be constrained) and consider a 20% 
uncertainty in this normalization, owing to 10% uncertainty 
in (is, ficiw and h each, over which we marginalize. We re- 
mark, however, that our results will be practically insensitive 
to the uncertainty on A. In Fig.(|9} we show the resulting 
constraints in the Q, m -w plane. Clearly, the Fksz method is 
more sensitive to w than 1Z, and their directions of degener- 
acy are also different, but both distinct to the directions of 
degeneracy corresponding to estimators based upon Large 
Scale Structure. Indeed, these estimators are restricted to 
the low redshift universe, and hence their sensitivity on w 
is very limited, giving rise to degeneracy dire ction almost 
parallel to the w axis, (see, e.g., Fig (11) in lEisenstein et alJ 
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FIG. 10. — Differential redshift contribution to the change in \ 2 when 
probing models with different w, approach a) is more sensitive to higher 
redshifts (thick line), method b) to lower redshifts (think line). 



noise level), then the resulting number of cluster pairs that 
can be formed would increase considerably, and this would 
reflect on an increased sensitivity to w . 



FIG. 9. — Contour plots showing the 1 and 2-cr joint confidence region 
(solid line) and the I— a marginalized (dotted line) in the Q m - w plane for 
the FK Z (method (a)) and 7Z 2 (method (b)) estimations. Method (a) applied to 
ACT (top) with CMB (WMAP 1st year) priors (dashed lines) yields a ~ 8% 
error on w. This error would drop below 5% if ACT covers 1,000 square 
degrees. Method (b) on a SPT-like survey gives a ~ 12% error on w after 
imposing the prior from WMAP's 1 st year data. 



(2005) or Fig. (13) in lTegmark etafl J2004)1. The constraint 
on w, marginalized over O m , is ~ 12% for F^sz and ~ 60% 
for 1Z. As the degeneracy direction in each case is different 
to that corresponding to CMB temperature measurements, 
a combination with analyses of CMB temperature data can 
break the degeneracy. When considering WMAP first year 
data, the marginalised error on w reduces to ~ 8%, ~ 12%, 
for Fksz and 1Z methods, respectively (thick dotted lines). 
As pointed above, these errors are dominated by the cosmic 
variance term present in ea.i'Ht. Therefore, by increasing the 
sky covered by the experiments, those errors should decrease 

aS OC 1/y/fsty. 

These results have been obtained after using a ST mass 
function, and assuming that the cluster density was uniform 
(i.e., we have have considered no biases). The error ampli- 
tudes and the orientation of the ellipses in Fig (|9j depend 
strongly on the number of cluster pairs that can be formed 
in each redshift bin, particularly at the high z end. If the 
actual sensitivity of future CMB experiments is such that the 
lower limit on detectable clusters can be relaxed (we believe 
the limit 2 x 10 14 h~ l M Q to be very conservative for ACT's 



We note that, a priori, the two methods are nicely com- 
plementary in their redshift ranges sensitive to w as shown 
in Fig. dlOl : we find that the F^ sz estimation (method a) 
thick histogram) is sensitive to w at high redshifts (z > 0.5), 
whereas 1Z for different dark energy models differs at low 
redshift and converges at z » 1 : the sensitivity of this method 
is localised mainly at low redshift (thin histogram). 

After parametrising the evolution of w as w(a) = wq + w a (l - 
a), with a = 1/(1 +z) the scale factor, we conduct a Fisher 
matrix analysis considering the parameter set [f2 m , wo, ficin, 
w a ] for the F^ sz method. The marginalization in the WQ-w a 
plane is given in Fig. ( II 1> . which shows this method should 
have some residual sensitivity on w a (a typical error of ~ 1 .75 
in w a ). Note that this error should be improved further after 
combining the two methods proposed here, provided that they 
are sensitive to different redshift ranges. 

7. CONCLUSIONS 

We have studied bulk flows in the large scale structure in 
the context of the kSZ effect and future high-sensitivity and 
high-resolution CMB experiments like ACT. Since the kSZ 
effect is only sensitive to the projected peculiar velocities of 
electron clouds, we have focussed our analysis on the angular 
correlation of radial peculiar velocities in galaxy clusters. 
We have provided an analytical expression for the angular 
correlation function of projected peculiar velocities in the 
linear regime, and interpreted it in the context of local and 
distant bulk flows. We have also presented an expression for 
the power spectrum of projected linear peculiar velocities. 

We have investigated in which redshift and for which 
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FIG. 1 1 . — Marginalisation in the WQ-w a plane for the ACT-like experiment. 
w a expresses the variation of w, which is parametrised as w(a) = wq + w a (l — 
a), with a the scale factor. Contours correspond to 1 & 2-cr joint confidence 
regions. 



cluster mass ranges large-scale bulk flows should be more 
easily detectable in future kSZ surveys, and computed the 
overall effect of the entire galaxy cluster population on 
the CMB sky. We have shown that the main contribution 
comes from Poissonian/random fluctuations of the number 
of clusters along the line of sight, especially at small angular 
scales. However, we find that the local bulk flow generates 
a signal which should be dominant at the very large angular 
scales of the quadrupole and octupole. Other terms associated 
to the coupling of velocity with density fluctuations give 
smaller contributions. We have calculated the kSZ signal for 
the cluster sample accessible by forthcoming experiments, 
and considered different sources of contamination which may 



limit our capacity to distinguish the kSZ signal from other 
components. 

We have presented two approaches to measure the kSZ 
signal and exploit its dependence on cosmological parameters 
such as the equation of state of dark energy (w). The first 
method is based on measurements of the kSZ flux and its 
redshift evolution. Its sensitivity to cosmology increases with 
redshift. For this reason it is suited for "deep and narrow" 
survey strategies. A data set such as that provided by ACT 
2 year observations with SALT follow up is well suited for 
an application of this method. ACT can detect the kSZ signal 
withaS/N~ 12atz~ 1. 

The second approach uses the ratio of kSZ and tSZ-induced 
temperature anisotropics. The cosmology-dependence is 
strongest at low redshift and hence this method is suitable for 
"wide and shallow" survey strategies. In this case we have 
considered an ACT-like experiment with larger sky coverage 
(1/10 of the sky) and increased instrumental noise. These 
are specifications similar to those of SPT with at least three 
frequency bands and combined with redshift determinations 
of the detected clusters. In this case the kSZ effect can be 
detected with a S/N of ~ 30. 

These methods can yield constraints on cosmological pa- 
rameters, in particular can constrain the equation of state for 
dark energy at the 10% level. The two methods are nicely 
complementary as they measure dark energy in different red- 
shift ranges, opening up the possibility to constrain dark en- 
ergy redshift evolution. 
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APPENDIX 
APPENDIX A 



Coupling the linear density and velocity fields 

In linear theory, the density contrast S(x) = (p(x)-p)/p is still much smaller than unity. This allows linearising the evolution 
equations and neglecting all non-linear orders, which makes the the evolution of each Fourier mode 5^ is independent from the 
other modes. 

In a Friedmann Robertson Walker (FRW) universe, the perturbed continuity equation reads 



dt 

which can be re-written in Fourier domain as 

d& 



^=-~ k ^ (A2) 
at a 

with a = 1/(1 +z) the scale factor. Throught this paper, the velocities are proper peculiar velocities. We must note that the 
expansion of the velocities in terms of its Fourier modes is given by 

/dk 
^3 v k exp-;kx, (A3) 
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and that, due to isotropy, the diferent components are independent. The statistical properties of each of them must be, however, 
identical. 

Coming back to ea.( IA2t . if we denote by T>g(z) the growth factor of the density perturbations, then we can express the compo- 
nent of vjj parallel to k (denoted here by Vjj) as 

^iHU-J-f. (A4) 

H stands for the Hubble function and z for redshift. 
The power spectrum for any component of the Fourier velocity mode hence reads 

P m { m ,dVs(z)\ 2 P m (k) 

Pvv{k)= [Hiz)—^) —. (A5) 

If we now combine 8^ and vq, we obtain 

<Vq> = (2nfS D (k-q)P dv (k) = (2nfS D (k-q) iH® * ^ PJ ® . (A6) 

Note that this is an imaginary quantity. In this work, we shall find this power spectrum either squared or in a convolution with 
itself, giving rise to negative power and indicating anticorrelation. 



Finally, we compute the power spectrum of the quantity 

#x,m) = (v(x) -m) = J2 S W ( ^ X K ! ) . 



(A7) 



where ni stands for the unitary vector connecting an observer with the position x, and the sum is over the three spatial components. 
Since a product in real space involves a convolution in Fourier space, the Fourier mode of reads: 



(A8) 



Having this in mind, the average product of two Fourier modes of when looking in two different directions ni and 112 is: 

dn 



*k,n,^n>(2^nk-q) 



cos 0i2 [Ps6*Pw] (£)+ (cos 12 +sin0 12 ) [Ps v *Pv8] (JO 



+ (27r)^"(k)(27r) j (5"(q) 



(2tt) 



3 PSv(u) 

(A9) 



012 = arccos(ni -112) is the angle separating ni and 112. Without introducing any loss of generality of ea. (IA9> . we have taken 
ni = (0,0, 1) and r\2 = (O,sin0]2,cos0i2). The symbol * denotes convolution and is present in the first two terms in brackets. The 
third term is non zero only when k = q = 0. 



Before ending this appendix, it is worth to make some remarks upon the redshift and mass dependence of the power spectra we 
have computed. The redshift dependence is explicit via the growth factors T>$ and T> v . Now, if we are interested in the peculiar 
velocity of a given cluster, then the (linear) velocity field must be averaged in a sphere of dimensions corresponding to the cluster 
mass. Something similar can be said about the density field: when studying the dependence of the number of haloes upon the 
environment density, the field 5 must be smoothed on scales which a priori are dependent on the mass of the haloes whose number 
density we are studying. In this scenario, all Fourier modes of the density and the velocity should be multiplied by the window 
functions corresponding to the scales within which we are averaging. This introduces a dependence on the masses of the clusters 
under study in P$s, P w an d Psv 



APPENDIX B 

The Angular correlation function of the projected velocities 

Let i and j be two components of the Fourier mode of the peculiar velocity field, so that 

i(v q )*\ = (2vr) 3 ,5 D (k-q) 6fj P w (k). (Bl) 
Using this, we can write the average product of projected velocities as: 

v(xi)-ni ) f v(x 2 )-n 2 
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k 2 dk 



P vv (k) W(kRi)W*(kR 2 ) / d<j)d(cos8)cos6 12 x 



(2tt) 3 

exp-/ k[x\cosO - jf2(sin 9 cos sin 6*12 + cos 6*12 cos 8)] . (B2) 

In this equation, the polar axis for the k integration has been taken along the direction given by ni. 6,<j) are the polar and 
azymuthal angles of k, and 612 is the angle separating the two directions of observations, ni and 112. If we now make use of the 
Rayleigh expansion of the plane wave, i.e., 

exp-ik x = (-0'(2/+ 1) ji(kx) P t (p), (B3) 

(where \i is the cosine of the angle between k and x and P\ are Legendre polynomials), and also the theorem of addition of 
Legendre functions (see, e.g., Gradshteyn & Ryzhik 8.794), then one ends up with 

v(xi)-m ] ( v(x 2 )-n 2 ) ) = V ^P- cos# 12 (- ft) [ k 2 dk P vv (k) W(kRi) W(kR 2 ) j t (k[ Xl -x 2 cosfl 12 ]) j,(kx 2 sm9 l2 ), 



even I 

(B4) 

where the factor T\ is given by 

(l-l)U IT 

Ti = —tr — ; — cos/ — , (B5) 
2'/2(//2)! 2' ' 

ji(x) are the spherical Bessel functions and the summation must take place only over even values of /. Note that in the limit of 
6*12 — * and X\ — * x 2 , this expression bec ome s eq. (0. 

At the same time, we can rewrite ea. JB2> as an expansion on Legendre polynomia. Indeed, if we denote by Xi and X2 the 
position vectors of the clusters, we can use the Rayleigh expansion of the plane wave for exp (/k • X! ) and exp (-/k • X2), and write 

v(xi)-m) (v(x 2 )-n 2 ^ = Y^(2l+\)(2l'+\){-it'' J j^PUk) cosfl 12 MkxOjdkx^Pip^Pri^), (B6) 

with /ix.y the cosine of the angle formed by vectors x and y. Note also that /ix, .x 2 = Mn, .n 2 = cos $12- Next we apply the addition 
theorem of Legendre functions on a spherical triangle formed by n 1 ,112 and k. As before, we take the polar axis of k to be aligned 
with Xi: 



v(x 1 )-n 1 j ^v(x 2 )-n 2 ^ =^(2/+l)(2/'+l)(-0'- / 'y pf n (i) cosfl 12 Mkx^j^kx^i^) 



2 P "'^k,x, ) p v (Mx, ,x 2 ) cos (m[0, - &]) 



(B7) 



Because dk = k 2 dk sin Bddd§\ , the integration in the azymuthal angle cancels the sum over m in the brackets. Finally, the prod- 
uct ^Xi,x 2 Pi(pxi,x 2 ) can be rewritten, via a Legendre recurrence relation, as a linear combination of Pi-\(pxi.x 2 ) an d Pm(MXi,x 2 )- 
After putting all this together, we find 

v(xi) • (y(x 2 ) ■n 2 )) = J2 C T' ^(cos0 12 ) (B8) 



47T 

I 



where the power spectrum multipoles CJ V are given by 



and the Br's are defined as 



CT' = ^j(/Bw + a+l)Bw) (B9) 



k 2 dk 

Bi = 4tt I -^Pwik) jiikx^jiikxi) (BIO) 

(27T)" 3 



16 



APPENDIX C 



The all sky kSZ correlation function 

Using ea.( !18> . we can write the average product of kSZ temperature anisotropies along two directions m, n 2 as 



kSZ 



T 



(nO^(n 2 )) = J dr 1 dndM 1 dM 2 dy l dy 2 f(Mi)f(M 2 ) W« fl! ( yi -r 1 )W^ s (y 2 -r 2 ) x 

v(y 2 ) 



{niMuyOi — -nj )n(M u y 2 ) 



(CI) 



Now we recall our model for the number of haloes of ea.(!16> to rewrite the ensemble average in eq. dCU as 

(n(M h y 1 )(^--n 1 )n(M 1 ,y 2 )(^--n 2 )) = n{M\ ,zi) 8 D {M\ -M 2 ) S D (j 1 -y 2 ) <jI,(M 1 ) + 



n(M u zi)n(M 2 ,z 2 ) 



v( yi ) 



11] 



C 

v(y 2 ) 



n 2 



dn(M u zi) 



OS 



dn(M 2 ,z 2 ) 



i=0 



OS 



(0(y 1; ni)0(y 2 ,n 2 )) 



(C2) 



5=0 



The first term is the Poisson term, and is zero unless both ni and n 2 are looking at the same cluster. The second term i s the 
velocity -vel ocity (vv) term, and the third contains the coupling of density and velocity studied in Appendix A. Plugging eq.(|C3) 
into eg. ( ICH and writing the integrands in terms of integrals in Fourier domain, one finds 



, 5T ksz 
" To 



(n0^(n 2 ))= / dndndMrdy, * ^ta-y.wq^y,, ^ (^V^MOnlM^,) + 

dr 1 dr 2 dM 1 dM 2 e- ikr ' +i ^n(M u zi)n(M 2l z 2 )Wr(Wrr f{M{)f(M 2 )t VR ' ^ 



dr\dr 2 dM\dM 2 



(2tt) 3 (2tt) 3 
dk dq -,kri+;qr 2 dn{M x ,zi) 



(2tt) 3 (2tt) 3 



<),S 



dn(M 2 ,z 2 ) 



,>=() 



ds 



Wr(W* as TT(Mi)T(M2) { 

=0 \ 



/ ^k,n, 



; q,n 2 



(C3) 



The integration on y l in the Poisson term generates a Dirac delta on k-q, whereas in the other two terms this Dirac delta arises 
naturally when one computes the ensemble average product of the Fourier modes of v and <fi. Hence, the integral on q dissappears 
and every integral ends up with a term of the form exp/k • (ri -r 2 ). We introduce now the Rayleigh expansion, yielding 



exp/k ■ (n -r 2 ) = £(2/ + 1)(2Z'+ 1)(-/) M ' ii{kn)j v {kr 2 ) Pi(^ kXl )Pr(p kr2 ), 



(C4) 



/./' 



and we make use (as in ea. JB7> in Appendix B) of the addition theorem of Legendre polynomia to express those Legendre 
polynomia having as argument /i^ r as a sum of Legendre functions of /i^ ri and /ir,,r,- Here, as before, we have aligned the 

polar axis of k along ri or ni . Also in this case, the integral of the azymuthal angle of k cancels the contribution of all Legendre 
functions having m ^ 0, so at the end one is left with 



,ST k 



kSZ 



To 



, JT ksZ/ ^2Z+1 [ f 

(ni )— (n 2 )) = ^— P/(cos0 12 ) ^ J 



(4 ' T) ~ lcdk dndr 2 dM x wr<w q n* f 2 (M,) Mkn)Mk r2 ) + 



(4tt) 2 k z dk 
(2tt) 3 



dr l dr 2 dM l dM 2 n(M u zi)n(M 2 ,z 2 )W g k as (W gas T f(Mi)T(M 2 ) 



P vv (k) 



ji(kn)ji(kr 2 ) cos6»i 2 + 



(4tt) 2 k 2 dk dn(M uZl ) 

5 — dr\dr 2 dM\dM 2 

(2?r) J 

(4tt) 2 k 2 dk 



(2tt) 3 



dr x dr 2 dM\dM 2 



86 

dn(M h zt) 



dn(M 2 ,z 2 ) 



OS 



8=0 96 
dfi(M 2l z 2 ) 

5=0 



rt=0 



W ga \W gas )* t(Mi)t(M 2 ) [Pdd*P w ](k) Mkn)Mkr 2 ) cos0 12 + 



<),s 



(>=() 



WftWFY r{M{)T{M2) [Pdv*Pr*](k) mn)ii{kr 2 ) ^sin0 I2 +cos0 

(C5) 

As in Appendix A, we have taken ni = (0,0,1) and n 2 = (0,sin#i 2 ,cos#i 2 ). Since different velocity components are not 
correlated, and we are only sensitive to the radial projection of the velocity, a cos#i 2 dependence appears in the vv term. For 
exactly the same reasons we obtained a cos$i 2 and a cos^ + sin^n dependence when computing the power spectrum of <f> in 
Appendix A. Note that, out of the three terms we found in that computation, we have dropped the last one because it is constant 
and introduces no anisotropy. Since the power spectra multipoles (C/'s) are projections on Legendre polynomia, such projection 
must be applied on cos6*i 2 and cosfl^ + sinf?^. The projection matrix for cos# 12 will be identical to that given by eq.( IB9> in 
Appendix B. For cos 9 l2 + sin 9\ 2 , it must be computed numerically. 
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Summarising, the power spectra corresponding to each of the terms considered here can be understood as a transformation of 
some vectors cf, (where X runs for Poisson, vv, dd-vv, dv-vd) by some linear applications Af v : 



cf = Y^Al v 4. 



(C6) 



For the Poisson term, A,f, sson is the identity, whereas for the vv and the dd-vv terms we find 



vv, dd-vv 



= 4tt 



151,,, a+i)$u' 



(2Z-1) 2 (2/ + 3) 2 

with Sfj is the Kronecker delta for i and /'. The projection matrix for the dv-vd term reads 



r+1 



The cf vectors for the Poisson, vv, dd-vv, and dv-vd terms are as follows: 



(C7) 



(C8) 



• Poisson: 



with Af(k,M) being 



• vv term: 



poisson = lj k 2 dkdM (Af(fc,M)) 2 , 

A?(k,M) = J drTavv[n(M,z)] 1/2 ji(kr)W k g 
cf = ~J k 2 dk (^A x ; v (k)j , 



(C9) 



(CIO) 



(Cll) 



with ( A^Ar)) 2 given by 



AJ v (k) \ = I dndr 2 dM l dM 2 f(M 1 ,zi)f(M 2 ,Z2) n(M u zi)n(M 2 ,z 2 ) W g k as (M l ,zi)W g k as (M 2 , z 2 )P n (M l ,M 2 ,zi, z 2 ,k) j,(kri)ji(.kr 2 ). 



dd-vv term: 



with 



Af~ vv (k) I = / dn dr 2 dM x dM 2 t(M u zi)t(M 2 ,z 2 ) 



„dd- 



-™ = 1 J k 2 dk^Af d - vr 



(C12) 
(CI 3) 



dn(M u zi) 



dS 



dn(M 2 ,z 2 ) 



<5=0 



dS 



Wg as (M u ziWr(M 2 ,z 2 ) x 



(5=0 



[Pdd *Pw](M u M 2 ,zi,z 2 ,k) ji (kri )j t (kr 2 ). 



(C14) 



dv-vd term: 



cf-vd = l j k 2 dk(Af v - vd 



(C15) 



with 



Af- Vd (k) = dr y dr 2 dM y dM 2 t(M u zi)t(M 2 ,z 2 ) 



dn(M u zi) 



OS 



dn(M 2 ,z 2 ) 



OS 



wr(M u zi)W g k a \M 2 ,z 2 ) x 



[Pdv * Pvd\(M x ,M 2 ,zi,Z 2 ,k)ji (kri )j, (kr 2 ) . 



(C16) 
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